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Explicit expressions are derived for the phase constant, the specific- 
group-delay constant, and the rms width of the impulse response for two- 
dimensional or square media having a transverse variation of index of 
refraction according to n = ni(l — \a u x u — |a r x r ), in which x is the 
transverse dimension, a u and a r are constants with \a r \ « a„, and 
(n — ni) <3C 1. Use is made of an approximation which the author has 
previously shown yields significant results. 

The results are applied to fibers with graded-index variation, clad by an 
additional medium of index n = n x (l — A). The ideal index gradient, a 
near-parabolic profile, gives delay distortion orders of magnitude less 
than for the conventional fiber with a step-change in index at the core- 
cladding boundary. However, it is shown that several forms of 5-percent 
error in the ideal gradient yield improvement of the order of 50 compared 
with the conventional clad fiber. The delay distortion is shown to be very 
sensitive to the exact index distribution in the vicinity of the ideal distribu- 
tion but increasingly insensitive to perturbations in the index distribution 
as that distribution departs more and more from the ideal. 

I. INTRODUCTION 

Optical fibers have assumed considerable importance for potential 
use as transmission media wherever wire pairs or coaxials are now used. 
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We give here some theory related to increasing the information 
capacity of such fibers. 

Conventional fibers have a core of uniform index of refraction n, 
surrounded by a cladding of slightly lower index of refraction. The 
cladding serves to isolate the outer fiber surface from the optical field, 
which is confined to the core, thus permitting the fiber to be handled 
and bundled into cables without affecting the information trans- 
mission. More recently, fibers having continuously graded index of 
refraction, with maximum on the fiber's axis and lower values at 
increasing radial distances, were proposed and realized in practice. 1-3 
Graded-index fibers have image-focusing properties 3 and provide 
modulated-carrier transmission with less delay distortion than con- 
ventional fibers having a step-index change at the core-cladding 
boundary. 4-6 Recent experimental work verified the potential of low 
delay distortion. 7 Current studies of graded-index fibers from Corning 
Glass Works show that total transmission losses under lOdB/kilometer 
and low delay distortion (approximately one nanosecond per kilometer) 
can be realized simultaneously in graded-index fibers. 8,9 This brings 
into focus the need for detailed knowledge about delay distortion 
related to the shape of the index gradient in the transverse plane. The 
following sections relate to that need. This study was done in parallel 
with that of D. Gloge and E. A. J. Marcatili. 10 The present paper 
presents an approximate analysis that may be extended to a wide 
variety of index distributions with closed-form solutions for the 
important wave-propagation constants. 

Section II summarizes the approximate method of calculation, 
yielding in turn the phase constant, the specific-delay constant, and 
the impulse response. 

Section III gives the solutions for the index distribution, 

n = ni(l - \a u x u — \a T x T ), (1) 

in which x is the radial coordinate, and a u and a r are arbitrary con- 
stants with a u > 0, a u y> \a,\, and 0.5 a u x u « 1. These conditions 
describe fibers of current interest. 

Sections IV through VI discuss cases of interest, taken as simplifi- 
cations of (1), which describe (i) an "ideal" index distribution, (it) 
variations around the ideal, and (in) other distributions which may 
result from convenient manufacturing processes. 

This paper is concerned with the delay difference between signals 
launched simultaneously in the various propagating modes ; the fiber- 
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output impulse response is derived assuming 

(i) all propagating modes are excited equally at the fiber input, 
(ii) the fiber structure is uniform along its length, yielding no 

mode conversion, 
(Hi) there are negligible losses or, equivalently, the same loss for 

all modes, and 
(iv) material dispersion due to variation of the bulk index of 

refraction versus wavelength is ignored. 

A very few modes very near cutoff are not accounted for; this is be- 
lieved to yield negligible error, and the same assumption was made by 
Gloge and Marcatili. 10 

It is found that the delay distortion of graded-index media is a very 
critical function of the index distribution in the vicinity of the near- 
parabolic distribution which gives the lowest delay distortion. For 
index distributions increasingly far from the ideal, the performance is 
less and less sensitive to changes. 

For fixed percentage error in fabrication of the index distribution, 
the delay distortion is linearly dependent on A, the fractional index 
difference between core center and the edge of the guiding region 
where the index becomes 

n = ni(l - A). (2) 

This is true even near the ideal index distribution, where previous 
work assuming no fabrication error indicated the delay distortion 
varied as A 2 . 11 

II. OUTLINE OF THEORY 

The approach taken here is based on Ref. 5, which gives an approxi- 
mate method for deriving the phase constant and other relevant 
quantities for wave propagation in the generalized media. 

We write the index n as a function of the transverse coordinate x, 

n = m[l + f(x)] (3) 

with |/(.r)| « 1. For guiding media, f(x) is predominantly negative; 
x = at the central axis. Here we assume f(x) is independent of the 
longitudinal coordinate. In accordance with Ref. 5 we derive a param- 
eter a e , which measures the transverse field width, using 

/(„.) =-o.i«5(^)'(£J- 
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Thus the phase constant /3 is given by, from Ref. 5, 

*-^M(£>+4 (5) 

The normal mode field may be considered composed of plane-wave 
components traveling at an angle a to the longitudinal axis; from 
Ref. 5, 

_ Xo (m + 1) 

a = -. (6) 

ni 4a, 

in which m is the mode number. The rath-order mode has (m + 1) 
extrema in the transverse cross section. The maximum angle that 
such plane-wave components can have is set by the fractional index 
difference A [defined in (2)3; any components with a larger than 

<w = V2A (7) 

exceed the critical angle for total internal reflection and are unguided. 
Thus eqs. (7) and (6) taken together establish a maximum modal 
index m max , with (m + l) ma x transverse extreme, controlled by a, 
and A. 

The specific group delay is 

<$ /m 

where <j = 2wf is the angular frequency. The range of values r can 
assume run from the value with m = to r max which is 

Tmax = t\ (m+l) mBX - (9) 

It is convenient to compare the specific group delay for the guided 
mode t to that for an infinite medium of index ni using 



'- ('-7| < 10) 

The work of Ref. 5 related specifically to two-dimensional wave- 
guides. We extend this here to three-dimensional guides in Cartesian 
coordinates. We note the relation between the propagation constants 
has the form 

cwSg = tf + & + z& , (11) 

where e r is the dielectric constant, j8 is the free-space phase constant, 
and /3 Z , j8»i, and /3 X 2 are the longitudinal and transverse wave numbers 
respectively. In our case, p x i and ^ l2 are small compared to fi t . The 
value of & depends only on the sum of the squares of the two transverse 
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wave numbers /& + $& ] further, for a square guide of width 2a e we 
can write 

p xl 2a e = (m a + 1)t, (12) 

X 22a e = (m b + 1)tt. (13) 

Thus the total modal designation is 



p = (m + 1) = ^(m a + 1) 2 + (m b + l) 2 (14) 

for the square guide. This replaces (m + 1) in the two-dimensional 
guide. The maximum (m + l) max can be reached with a variety of 
field distributions given by various values of m a and m b in (14). 

To compute the impulse response we first note in solutions for 
specific group delay r from (8) that all combinations of modes having 
the same m or p in (14) have the same specific group delay. Hence 
the fiber output at a given delay associated with p will be P{p), for 
equal power into the fiber in each mode, where 

P( P )=^d P (15) 

and p is given by (14). The fiber impulse response as a function of 

time t is 

dt 



Output = P(t) = ^ = ^ 



dp 



(16) 



The relation between t and p is obtained from (10) and (8). The range 
of t over which (16) is valid is found by inserting the minimum and 
maximum values that p may assume, 

p mi „ = V2, (17) 

p ma x = (m + l) max . (18) 

For many purposes the value of t miD corresponding to p in in may be 
taken as zero since p ma x ^> Pmin ; thus 

* mi „^0, (19) 



t max = (t^-^) 



(20) 



The effect on transmission system performance is approximately the 
same for various shapes of the impulse response if the rms width is 
the same ; 12 this applies in the region where the fiber impulse response 
degrades the system signal-power requirements by only 1 or 2 dB. We 
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find the desired second moment using 

A = f tm " P(t)dt, (21) 

Jo 

1 t 'max 

T = j j o tP(t)dt, (22) 

1 /"'max 

<r* = j J t*P(t)dt - T\ (23) 

III. A GENERAL SOLUTION 

We now give the results for the index distribution according to 
eq. (1). The guide has width 2a. We specify that at x = a, f(x) = — A, 
leading to eq. (2) and 

«, = ^, (25) 

in which 5 may be either positive or negative but |5| <3C 1. Using eq. 
(4), we find a e , 



1_ 

«2 



n\a u 1 2 "«+ 2 > 2a r n 2 

H r "T" TTTTTT^ ' (26) 



°« ( (« + 2X..3036-XS I *■ 1 " < " + " 



0.303& 4 Ag 
This leads immediately to /3, 

2wni tt(w + l) 2 



£ = 



16(0.1515) 2 /<»+ 2 >^ m ^ 5 2 - 5 V 



[(1 - 5)A] 2 '<"+ 2 > / \ y«-2)/(«+2) 



- 2.592 X 



fl 2u/(u+2) 

(0.1515) r/(u+2) (m + l)2A<«+2-'-)/(»+2)5 



/ W -4- 2 ft \ (2u+4-2r)/(u+2) 
+ 2) (^5^j (1 - *)'"«+»*' 



/(u+2) 



*(£) ' (27) 

and to the number of propagating modes N, 

N = - (0.7757) 2 /"(ni/ca) 2 A, (28) 

7T 

in which k = 2ir/\<,. The maximum modal index has two useful forms, 
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derived from (6) and (7), 

(m ( +2 5?%+*> = 32 (°' 02424 ) 2/(u+2) ( ™ J" h2> A" /(u+2) , (29) 

ro max ^5.657(0.7757) 1 '«(^) A*. (30) 

The specific group delay from (8) is 

_ m I (u-2) _ /(u+2) (m + l) 2 (m + 2.5)Vg+» 

T " 7 J 1 + fiT+2) AU j (m + l)Lx (w + 2.5)«'<«+» 

g (m + lY (m + 2.5)^ (n+2) 1 

+ ^ A (1 _ 5)r/(.-H) ( ?n + 1)52 (m + 2.5)< 2 «+«-">'<» +2 > [ ' ^ 1; 
where Q is 

Q = 2.578(0.7757)'-/" ^j^ - (32) 



We can simplify r for m » 1 to 

i ■ • («zJ) A( i - a) ./(.« >/^Y u/(u+2) 

(tt + 2) V ' V^max/ 

+2) 



r = 

C 



5 / m V r/(u - 

Using (21), (22), and (23) we find the impulse response is 
P( A I ,2/« I ir(n + 2) «"">'"*&., 1 

/ W l«^-« | 4m ( w _ 2 )(«+2)/ uA (u+2)/ U ( 1 _ 5 )2/u I 



(33) 



(34) 



4w 2 (n - 2)< 2+r+u)/u (l - 5)('+2>/»A (r+2)/n 
and the rms width of the impulse response is 

n * M on/ 1 (m_+2) V rtW 

, - 7 A(u - 2) {^-^ - (2w + 2)2 J (35) 

for the case where a r = 0. The minimum allowed value of t in (34) is 

niA (m - 2) 2 (m + 2.5) 4 m / ft < x " +2) f . 

' min ~ c (tt + 2) (2.5) 4 '<« +2 > (m + l)Lx 

Equation (29) can be used to eliminate the last factor in (36). For 
the first term of (34), representing the major response due to a u x u in 
(1), the maximum value of t is 

_ niA (w - 2) r „-v 

«„.«(«)- — (l4 + 2 )" ^'' 
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For the second term of (34), representing the perturbing term a^x, r 
in (1), the maximum value of t is 

tmax(r) c Q _ g\r/(u+2) V°°/ 

Equation (34) is not valid for u = 2; for the later case the impulse 
response is 

4096(1 - 5)"«ro5& I*- 01 '. (39) 



P(0I— i = -i (4 - r)/r - 



32 r / 2 (41.254) w-')/»(2r - 4)5A 



IV. THE NEAR-PARABOLIC INDEX DISTRIBUTION 

Letting u = 2 in the equations of the preceding section yields the 
near-parabolic index distribution. As shown by (35), the impulse 
response has zero width when a r = in the approximation made here. 
In the cylindrical fiber there is no index distribution which gives zero 
delay distortion among all the various modes. There is a distribution 
of index which minimizes the delay distortion, and we now evaluate 
this condition. We can use the above theory to evaluate the cylindrical 
waveguide by noting the two limiting conditions already known for 
low dispersion. Take the form 

n = ni (i - i a 2#2 + & 4a 4#4 + . . . ) (40) 

in which R is the radial transverse coordinate. In two earlier papers 8,13 
it has been pointed out that the value of 6 4 must be different to produce 
no dispersion for meridional rays versus skew rays; the difference in 
6 4 was found to be 5. Kawakami and Nishizawa 6 found that 64 must 
be \ to give no dispersion for skew rays and must be \ to give no dis- 
persion for meridional rays. We can visualize minimizing the dis- 
persion for one ray type, and thus experiencing a maximum dispersion 
corresponding to a change in 64 of \. We do this in the approximation 
used here by setting u = 2, r = 4, and making 

a 4 =-ia| (41) 

corresponding to no dispersion for meridional rays. From (24), (25), 
(31), and (41) we find for this "ideal" index distribution in the round 
fiber the specific group delay 

(m + l) 2 , 



= V± 

c 

and 



1 - 0.26 



(42) 



(m + 1)L« 
* max =-0.26^A>. (43) 
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This agrees reasonably well with the value — riiA 2 /8c arrived at by 
Gloge and Marcatili 10 using an entirely different analysis for an 
optimum index distribution denned differently.' The rms width of 
the impulse response corresponding to (41) and (42) is 

a = 0.752 ^ A 2 . (44) 

c 

In contrast with this, the simple step-index fiber [represented by 
a r = and n — >°o in (1)] has an rms impulse response width of 

L^A. (45) 

12 c 



Thus the "ideal," corresponding to (42) and (43), is smaller by a 
factor of 0.26A. Since A ~ 0.01 the ideal graded-index fiber has an 
impulse response narrower by a factor of about 400 than that for the 
conventional fiber. 

The fourth-order term represented by (41) corresponds to 

5=-§A. (46) 

This is very little different from the simple parabola — not enough to 
see on Fig. 1 where curve ® represents both 5 = and (46) for 
A ~ 0.01. More importantly, (46) and (41) imply that the fourth- 
order term decreases in size relative to the second-order term as A 
decreases. The inaccuracies of material processing are likely to prevent 
this as A becomes small. A more probable limit is a fixed value of 5 
in (25) as A changes. This results in a specific group delay 

' 0.3895 (m+ l) 2 A \ (47) 

and an rms width for the impulse response 



c 



, = 0.112^A 7 - r ^- v (48) 

c (1-5) 



* If instead of minimizing the delay distortion for either the skew rays or meridional 
rays we had minimized the delay distortion for an index equation (40) at the mean 
of the index values giving minimvun distortion for the skew and meridional rays — i.e., 
at b = 5/12— then the maximum departure in b for any mode corresponds to a 
change in b =±1/12. This corresponds to a 4 =±1/Ga\, leading to 

r=g{l±0.13 / W + y A- 

cl (m + iy mmM 

in place of (42). The coefficient 0.13 corresponds almost exactly to the result of 
Gloge and Marcatili, 10 but the present analysis indicates twice the total delay spread 
predicted by Gloge and Marcatili due to the ± sign. 
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Fig. 1 — Normalized index of refraction versus transverse coordinate (x/a) for the 
following parameters in eqs. (1) and (25) : curve ®, u = 2, & = 0; curve ®, u = 2, 
r = 4, 5 = 0.05 ; curve ®, u = 4, 5 = ; curve ©, w = 4, r = 2, 5 = — 1 ; curve ®, 
u = 8, 5 = ; curve ®, u = <*> , 6 = 0. 



For 5 = 0.05, a becomes 0.00591 niA/c which is narrower than for the 
step-index fiber by a factor of about 50 independent of A. 

We note from (39) that the impulse response for the ideal index 
distribution perturbed by a fourth-order term (r = 4) is a rectangular 
pulse, shown as curve © in Fig. 2. However, if the perturbation were 
sixth order, r = 6, the impulse response would vary as t~*. Other values 
of r give other impulse-response shapes, which we discuss further in 
the next section. 

Finally we note from (47) that the impulse response due to the 
fourth-order perturbation of the ideal distribution may either lead or 
lag the impulse at r = ni/c, depending on whether 5 is positive or 
negative [see (1) and (25)]. Similar effects due to the sign of 5 are 
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Fig. 2 — Impulse response P(t) versus t/t^^ where t = t - n\/c and t is given by 
(31) and (33). The conditions are the same as defined under Fig. 1. 



found for perturbations of the nonideal index distributions and may 
be seen in (31). 

In Section VI there is a discussion of Fig. 8 which shows the effects 
of perturbing the parabolic index distribution in several ways. 

V. DISCUSSION OF THE INDEX DISTRIBUTION, n = n,(1- 1 /ia„x") 

In this section we discuss the distributions obtained by setting 
a r = in (1) which mean 5 = in (24) and (25). 

The total spread in specific group delay for all modes t maI , u is given 
by (37), which gives a null value when u = 2. As already discussed, 
this is a simplification in the vicinity of the "ideal" distribution. 
However, (37) gives a valid representation as u departs significantly 
from the value 2. Figure 3 shows the variation in t m ^, u versus u. 
The behavior in the vicinity of u = 2 is a form of singularity. For 
5-percent error in u from the ideal, « m8X , u ~ 0.025niA/c. This com- 
pares with 0.02447uA/c for the modal delay spread from (47) for 
5-percent (at x = a) fourth-order perturbation of the ideal. We con- 
clude that it is not particularly important how the ideal is perturbed. 



GRADED INDEX FIBERS 187 




ORDINATE >0 

ORDINATE <0 



1.5 2.0 2.5 3.0 3.5 4.0 4.5 

u 

Fig. 3 — Normalized modal group-delay spread versus u for a T = in eq. (1). 



The value of t max , u from (37) and plotted in Fig. 3 is not very 
sensitive to the value of u at values removed from u = 2. The reduction 
in rms width of the impulse response is illustrated in Fig. 4. For u = w 
(the conventional step-index fiber) the value of <r/(niA/c) is 1/Vl2 or 
0.289. This is reduced by a factor of 2 for u near 6, and by a factor 
of 4 for u near 3.5. These results are identical to those of Gloge and 
Marcatili. 10 

The shape of the impulse response, given by (39) for u = 2 and 
by (34) for u away from 2, is plotted in Fig. 2 for several cases of 
interest. 

The number of modes which can propagate, given by (28) for the 
square fiber, is plotted as a function of u~ l in Fig. 5. For comparison 
the corresponding quantity for the round fiber from Ref. 10 is also 
plotted. The ratio is 4/7r at u = °° , and near 2 for u = 2. The approxi- 
mations made here are seen to be good, though not perfect. 

VI. PERTURBATIONS OF THE GENERAL INDEX DISTRIBUTION 

We discuss now some of the results for the perturbed index distribu- 
tions, eq. (1). We recall the solutions have been obtained assuming 
I Or | « a u or \h\ « 1. 

The solutions for the specific group delay, given in (31) and (33), 
contain the quantity Q as a factor in the perturbing term. The factor 
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Fig. 4 — Normalized rms width of the impulse response versus u for Or = in eq. (1). 




Fig. 5 — Total number of propagating modes versus u l for a r = in eq. (1). 
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Q, given by (32), contains the principal effect of the exponents u and 
r on specific group delay. Figure 6 shows how Q varies with r for the 
special case u = 2. The region verj r near r = 2 is in question since we 
know the "ideal" index distribution differs slightly from u = 2. Else- 
where, the results should be significant and may be used in (34), 
(31), and (33). 

More general curves for Q are given in Fig. 7. We observe that 
when r > the maximum value of Q is not very dependent on u but 
the most sensitive region of r (giving largest Q) does depend somewhat 
on u. An intuitive feel for the changes in the index distribution which 
correspond to some of the curves in Fig. 7 can be obtained by examina- 
tion of Fig. 8. Figure 8 shows the normalized index n versus transverse 
coordinate (x/a). The curve labeled r = 2 corresponds to the pure 
parabolic distribution. The other curves correspond to 5 = 0.05 with 
various values of r in the perturbing term and u always equal to 2. 

We note in Fig. 7 that, at u = 2, the value of Q at r = 10 is much 
larger than at r = °° . This may seem surprising, since a step-index 
change occurs at (x/a) = 1 .0 when r = °° . Below (x/a) = 1 the 
r = oo curve in Fig. 8 corresponds to a pure parabolic gradient between 
the ordinate equal zero and — 0.95A. 
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Fig. 6 — Q versus r for the near-parabolic index distribution. 
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40 60 80100 



Fig. 7 — Q versus | r | with u as a parameter. 

We also note in Figs. 8 and 7 that dips in the index distribution 
near (x/a) equal zero (curves for r = -0.4 and -0.1) yield values of 
| Q | comparable to those for r in the range 4 to 20. 

VII. CONCLUSION 

The above analysis provides an approximate solution for the delay 
distortion to be expected in a wide variety of graded-index fibers, 
representable by (1) with \a r \ « a u . 

In general, gradual tapering of the index between the center of the 
fiber and the outer support provides reduced delay distortion. Only 
in the vicinity of the near-parabolic distribution is the performance 
highly sensitive to the exact index distribution. The "ideal" near- 
parabolic distribution provides a potential reduction in delay distortion 
of several hundred times compared to the step-index distribution of 
the conventional clad fiber. With an accuracy of the order of 5 percent 
in achieving the "ideal" distribution, the reduction in delay distortion 
is on the order of 50. 
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-r = 10 OR to 




Fig. 8 — Normalized index versus transverse coordinate (x/a) for several per- 
turbations of the parabolic index distribution. 
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